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SYMBOLS 


a local speed of sound, m/sec 

CD drag coefficient of a profile 

CL lift coefficient of a profile 

CM moment coefficient of a profile (with respect to quarter-chord point) 

CP pressure coefficient 

M local Mach number 

P pressure 

R rotor radius, m 

r/R relative radius of blade station 

t time, sec 

wind velocity in the wind tunnel, m/sec 
x/c chordwise relative position on a profile 

ot Q incidence, deg 

y ratio of specific heats (1.4 for air) 

y advance ratio (U/ftR) 

p density 

<j> velocity-potential in the inertial frame 

[ <J> ] jump in velocity-potential 

^ blade azimuthal position, deg 

Q angular velocity of the rotor, rad/sec 
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SUMMARY 


A numerical method is presented for calculating the three-dimensional unsteady, 
transonic flow past a helicopter rotor blade of arbitrary geometry. The method 
solves the full-potential equations in a blade-fixed frame of reference by a time- 
marching implicit scheme. At the far-field, a set of first-order radiation condi- 
tions is imposed, thus minimizing the reflection of outgoing wavelets from computa- 
tional boundaries. Computed results are presented to highlight radial flow effects 
in three dimensions, to compare surface pressure distributions to quasi-steady pre- 
dictions, and to predict the flow field on a swept-tip blade. The results agree well 
with experimental data for both straight- and swept-tip blade geometries. 


iNl'KOuUCi lON 


The exact flow fields around helicopter rotor blades are generally acknowledged 
to be complex, unsteady, and three-dimensional. Problems arise because of strong 
flow separation on the retreating blade and because of transonic effects on the 
advancing blade. The present study is concerned with the advancing blade, the tip of 
which may enter the transonic flow regime and the aerodynamic coefficients of which 
may change appreciably as a result of transient shock-wave movements. Accurate pre- 
dictions of these shock motions are essential preconditions to further improvement of 
the helicopter rotor blade. The difficulty and expense of obtaining such detailed 
flow information experimentally are sufficient reasons for developing a computational 
method for characterizing the flow. 

All previous computer codes for handling these flows have involved small distur- 
bance on mean surface approximations (refs. 1-5). The most advanced of these codes 
is an alternating-direction, implicit solution scheme for the low-frequency, unsteady, 
transonic small-disturbance equations. In the present study, the flow physics is 
modeled by the unsteady, transonic full-potential equations; this modeling allows 
more realistic specification of the blade geometry, and should provide better reso- 
lution of the associated flow phenomena. The numerical method for two-dimensional 
airfoils, developed in reference 6, was extended to calculate three-dimensional, 
unsteady rotor flows. 

This is Part II of a series of planned publications under the same general title 
"Transonic Flow Analysis for Rotors." 


HELICOPTER ROTOR FLOW FIELD 


Before discussing the details of the numerical model, it is worthwhile to dis- 
cuss, in general terms, helicopter rotor flow fields. The helicopter rotor blade 
generates an extremely complex flow. The combination of forward speed and uniform 
rotation causes each blade section to encounter a sinusoidally varying free-stream 
velocity. In addition, to avoid rolling moments, blade incidence is varied nearly 
sinusoidally throughout each blade revolution. The result is a very complicated 



aerodynamic environment. Figure 1, taken from reference 7, was obtained by a calcu- 
lation based on an acceleration potential method; it shows contours of constant Mach 
number, constant-sweep angle, and constant incidence over the rotor disk. 


GOVERNING EQUATIONS 


In order to develop an accurate, efficient, and reliable computer code that can 
be used for engineering purposes, the flow modeling is based on potential-flow theory. 
The shock wave is therefore assumed to be weak enough so that the entropy jump across 
it is negligible. In other words, the flow is assumed to be inviscid and isentropic 
over the entire flow field. A velocity potential <J> exists for the flow and is 
described in a frame of reference that is at rest relative to the undisturbed air. 

In this inertial frame, the complete equation for the velocity potential is 

<P tt + [(V*) 2 ] t + V<t> • v T-i- <fc).l = a 2 V 2 (J> 


where a is the local speed of sound. Bernoulli T s equation, relating a and 



(V *) 2 + 




is 


where a^ is the sound speed in the undisturbed air, and y is the specific heat 
ratio (1.4 for air). 


If the blade geometry and location are described by S(r;t) = 0, where r is 
the position vector in the inertial frame, then the boundary condition at the blade 
surface is 


S + % • VS = 0 

For further analysis it is more convenient to implement this surface boundary condi- 
tion in a moving frame of reference in which the blade location is fixed. Let primed 
variables refer to the inertial frame F ! and unprimed variables refer to the blade- 
fixed moving frame F. Suppose that at time t, the two frames are coincident and 

that F is moving relative to F' with a linear velocity U and an angular velocity 

Then, at that time t, the position vector r of a particular fluid particle is 
the same for both frames. If a point P is rigidly fixed in F ! , it appears to an 
observer in F to move with velocity ^ = - ($ 4- fi x r) . Thus, the velocity of fluid 

particle at P in F is q = ^ + V. The rate of change of $ at P is measured 

by an observer in F as 


. f = * + V • V(f> 


T t ' T t 

The potential equation in the moving frame F is given by 




<P tt + 2V • 7<f> + (V • V) (V • V<j>) + V<j) • V fc + 2V<p • V<j> t 


+ (^ • ^) [ (V<j)) 2 ] + V(f) 


V jl (V,) 2 ] = 


a 2 V 2 (J> 
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and Bernoulli’s equation is 


i 2 a 

->■ 1 o a 00 

1> t + V • V* + 2 (V<f>) 2 + y~r~[ = T^TT 

Let the moving frame F be described in a Cartesian coordinates system in which 
x, y, and z represent the chordwise, vertical, and spanwise directions of the blade, 
respectively, and the origin of which is at the center of rotation. In the inertial 
frame, let the advance velocity I? lie in the (x’,y T ) plane and let it make an 
inclination angle a Q with the negative x’-axis direction, and let the angular 
velocity fi be in the positive y’-axis direction. The velocity ^ owing to the 
motion of the moving frame F has the components 

V = ttz + U cos a sin 

1 o r 

V 0 = U sin a 

2 o 

and 

V 0 = -fix + U cos a cos 

3 o 

where \jj is the azimuthal angle of the blade .(^ = 180° for forward flight direction 
in the inertial frame) . 


The potential equation in Cartesian coordinates is 

d> 4* 2q 6 + 2q d) + 2q (f) 

y tt ^l y xt H 2 Y yt M 3 y zt 

= • (a 2 - q 2 )(f> 4- (a 2 - q 2 )<|) + (a 2 - q 2 )<j> - 2q q 6 - 2q q d> 

Y XX M 2 /Y yy v Y ZZ n l n 2 Y xy 13 

+ (ft 2 x - 2W cos a cos + (ft 2 z + 2ftU cos a sin ^)<() 


xz 


2q q d> 
H 2 M 3 Y yz 


( 1 ) 


where q x , q 2 > and q 3 are the velocity components of local fluid particle in the 
moving frame and are specified as 


q = <j) 4 V 

Y x l 

q = d> 4 V 

4 2 Y y 2 

and 

q 3 = + V 3 


Bernoulli’s equation in Cartesian coordinates is 


<(k 4 V A 4 V A +Vi 
Y t l Y x 2 Y y 3 z 

Several boundary conditions are 
In the near-field, the flow tangency 


+ 4 (<f> 2 + (j> 2 4- <j> 2 ) 4- 

2 VT X Y y Y Z 


Y - 1 y - 1 


( 2 ) 


necessary to complete the boundary value problem, 
to the blade is described by the expression 


q • n = 0 
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where n is the unit vector normal to the blade surface. A simple rectilinear wake 
model is used here. The wake that is shed from the trailing edge is assumed to be a 
vortex sheet that is a smooth continuation of the trailing edge. Across this vortex 
sheet, normal and binormal velocity components and pressure are assumed to be contin- 
uous. The jump in potential determined at the trailing edge of each spanwise profile 
is then assumed to propagate to infinity according to the wake condition 

[4> t ] + q m • = o O) 

where q m stands for the mean of upper and lower velocities at each point on the 
vortex sheet. At the far-field, the boundary condition can be formulated as a 
Neumann condition so that the time-derivatives of the velocity potential along the 
characteristics of the flow equation in each spanwise plane vanish. Mathematically, 
it is subject to the condition 

<f> t + = 0 (4) 

where q a is the outgoing eigen-velocity of the flow equation in each spanwise plane 
(ref. 6) . 


MESH SYSTEM 


A parabolic sheared mesh system (refs. 8-12) is used in the present analysis. 

The blade surface is conformal with a coordinate surface of this system. Thus, it is 
easy to include the blade surface condition in a finite-difference calculation. The 
mesh system is generated by a series of transformations from the physical space to 
the computational domain; it is briefly described here in the interest of 
completeness . 

The shearing transformation- The shearing transformation, 

X = x - x g (z) 

y = y - y s (z) 


and 


z = z 

shears out the blade sweep and dihedral. The point x s (z), y s (z) is the center of 
the circle passing through three points near the leading edge of the profile at each 
spanwise station. 

The scaling transformation- The scaling transformation, 

x = x/SCAL 
y = y/SCAL 

and 

z = z/SCALZ 
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accounts for the scaling between the physical space and the computational domain; 

SCAL and SCALZ are real variables, which accounts for the scaling between the physi- 
cal space and the computational domain. 

The square-root transformation- The square-root transformation, 

(X x + iY,) 2 = 2(x + iy) 

and 


maps the entire blade surface to a shallow bump Y 1 = S(X l9 Z 1 ) near the plane 
Y x = 0. 


The second shearing transformation- The second shearing transformation, 

X = X x 

Y = Yn - S(X x ,Z x ) 

and 

Z = Z x 

reduces the blade surface to a portion of the plane Y = 0. 

The stretching transformations- The stretching transformations are introduced to 
render the computational domain finite. For example, 

Y = — , 0 < a < 1 and b > 0 

(1 - Y 2 ) a 

is used to map the planes Y = ±°° to Y = ±1. Similar transformations are used out- 
board of the blade tips in the Z-direction and downstream of the trailing edge in 
the X-direction. 

The flow equation (1) in the computational domain is very complicated and will 
not be repeated here (the reader is referred to ref. 11 for details). There is a 
degenerate case. For points on the continuation of the singular line outboard of the 
blade tips, where the Jacobian vanishes, the potential equation reduces to a simple, 
two-space-dimensional Laplace equation. 


FINITE-DIFFERENCE SCHEME 


The time-marching scheme developed in reference 6 for two-space-dimensional 
problems will be extended to three-space-dimensional problems. The scheme is of 
three time-levels, with time-levels n - 1 and n predicting time level n + 1 . 

There is no iteration in any of the time-marching procedures. Basically, the velocity 
components at each grid point are evaluated by using central differencing at time- 
level n. The time-derivative term in Bernoulli’s equation is evaluated by backward 
difference. 
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It is convenient to have the following difference operators: D central differ- 

ence operator, 5 upwind difference operator, and D-type difference operator. To be 
specific, the definitions of operators acting in the x-direction are given by the 
following: 

1. Central difference operator: 

D xh = ( h+i - h-i )/(2Ax) 


2. 


3. 


Upwind difference operator (forward or backward): 


q D <j> . = \ 
n i x T i 1 


~ if 

[_ q i * ^i+i ” t f> i )/( Ax ) if 


Type-dependent difference operator : 


for subsonic points 


q x > o 
q x < o 


D 

X 


X 


for supersonic points 


The second-order difference operators are defined as combinations of the above dif- 
ference operators; namely, D X y = D X Dy , etc. 


Flow Equation 

The potential-flow equation in canonical form is 

<j> + 2 (q • V) <j> = (a 2 - q 2 )c + a 2 (V 2 cj> - $ ) + first-order terms (5) 

u t t s s s s 

where s is the local streamwise direction unit vector. Its finite-difference 
approximation can be described as the following. For the left-hand side of the above 
equation, upwind differences are used for all spacial derivatives, and central dif- 
ference is used for all temporal derivatives; namely, 

*tt + 2q At + 2q 2 V + 2q 3+zt 

S (4> n+1 - 2<p n + 4» n-1 )/(At) 2 + (q^ + q 2 D y + q^) (<j, n+1 - <f, n-1 )/(At) 

Here, we add some convection viscosity through upwind differences which stabilizes 
the scheme but introduces some dissipation. The desired ability to capture shocks 
automatically requires the introduction of viscosity in the differencing of the right- 
hand side of equation (5). This is done by means of Jameson’s rotated difference 
scheme (refs. 8-11). That is, for contributions to <j> ss , type-dependent differences 
are used for all spacial derivatives, and the term <j> n is replaced by the arithmetic 
mean of <f) n+1 and (f) n_1 : 
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+SS ■ r? <<\ x + <+,y + q 3+zz + 2q iV*y + 2 V3+ XZ + 2q 2 < >3*yz ) 


T (q z D + q 2 D + q*D + 2q q D + 2q q D + 2q q D )(<j>“’* + *) 

2 XX 2 yy ^3 ZZ n X n 2 xy XZ M 2^3 yz Y 


n+l 


^n-l. 


2q 2 


yy 


xy 


For contributions to V 2 (j> - <j> ss , central differences are used for all spacial 
derivatives, and the term (f) n is again replaced by the mean of <f) n+1 and (J) 11 " 1 . For 
the first-order terms in the right-hand side of equation (5), central differences are 
used for all spacial derivatives • Finally, the finite-difference approximation to 
the whole flow equation can be written as 

(* n+1 - 2<p n + tp n_1 ) / (At) 2 + (q x D x + q 2 D y + q^) (<t> n+1 - ^"^/(At) 


= (a - q )/(2q )(q D + q 2 D + q D + 2q q D + 2q q D + 2q a D ) 

\ M XX M 2 yy ZZ 2 xy 1^3 XZ M 2 M 3 yz 


• (<()“' A + <f> n b + a z /(2q z )[(q 2 + q 2 )D xx + (q 2 + q 2 )D + (q 2 + q 2 )D 


ZZ 


- 2q q D -2qqD -2qqD ] (<(. n+1 + <j> n 1 ) 
n l n 2 xy 13 XZ n 2 n 3 yz Y 

+ (ft 2 x - 2^2U cos cos i|j)D x ({> n + (Q z z + 2Q\J cos a Q sin ijODyJ) 11 

The discretization error associated with the finite-difference approximation is 
of second order in space for subsonic points and first order in space for supersonic 
points and second order in time. The system of algebraic equations generated by the 
finite-difference equation is large. However, this equation can be factored within 
the same order of accuracy. The following factorization has been tested and found to 
be stable with time-steps much larger than the time-step allowed by the CFL condition 
for an explicit scheme. 

Let M = q/a; then 


L = 1 + Atq D + 


X 


1 X 


At q 2 D - (q 2 D + q 2 D + q 2 D ) 

'll -,2 xx n 2 XX n 3 XX 7 


q xx 


2M 


L = 1 + Atq D + — “ — q 2 D 
y n 2 y ° ° ■ 


and 


L z - 


on- - ~~r- (q 2 D + q 2 D + q 2 D ) 

2 2 yy 2 m 2 1 yy 2 yy 3 yy 


1 + Atq D + q 2 D - (q 2 D + q 2 D + q 2 D ) 
M 3 z 2 ^3 ZZ 2M 2 1 2Z 2 zz 3 zz 


Then the approximate factorization of the finite-difference equation can be written as 


LLL <}> n+1 = RHS 
x y z 


Here RHS consists of the terms involving <|> n and <f> n in the above finite- 
difference equation. This factorization reduces the large complicated matrix inver- 
sion problem to three one-dimensional problems. The solution algorithm can be 
expressed as 
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Each of the above steps required a penta-diagonal matrix solver because second-order 
upwind difference operators are used in each of the difference operators L. 

Radiation Boundary Condition 

The artificial radiation boundary condition (eq. (4)) is of the form 

<b + ud> + vd> + wd) — 0 
Y t Y x T y T z 

which is approximated by 

($ n+1 - <f> n_ 1 ) / (2At) + (uD x + vD y + w$ 2 )(<j) n+1 + <t> n_1 )/2 = 0 

Again, it can be factored as 

(1 + Atuf> x ) (1 + Atvf) y ) (1 + AtwD z )(J. n+1 = RHS 

The algorithm consists of the following three steps: 

(1 + AtuD x )X = RHS 

(1 + AtvD )Y = X 

y 

and 

(1 + AtwD z H n+1 = Y 
Wake Condition 

If the wake of each blade section is assumed to lie on the x-axis, the wake 
condition is 

[<)>.] + u [<(>]= 0 
LY t m Y x 

where u m is the mean of upper and lower wake velocities. The finite-difference 
approximation is given by 



Let b = i^At/Ax; then 



1 + b 


Hence, once the jump at the trailing edge has been found, the jump in the entire wake 
can be calculated from the above formula. 


RESULTS AND DISCUSSION 


A new computer code that contains both quasi-steady and unsteady modes has been 
developed. The quasi-steady mode is based on the theory reported in reference 11, 
and the unsteady mode is based on the above theory. Good results can be obtained on 
mesh containing 128 x 24 x 32 grid points in the chordwise, vertical, and spanwise 
directions, respectively. On the Ames 1 Cray 1-S computer, it takes about 40 sec of 
CPU time for the quasi-steady mode to calculate each azimuthal position and about 
30 min CPU time for the unsteady mode to calculate a half rotor revolution, from 0° 
to 180° (with 0.25° per time-step). Further improvement in computing speed can be 
realized through a more complete vectorization of the code. The capability of the 
computer code is demonstrated through the following studies. 


Radial Flow Effect 

First the quasi-steady mode was used to obtain radial flow effect. Two quasi- 
steady calculations were performed on the same ONERA, nonlifting, straight-tip rotor 
blade (ref. 12) at azimuthal angles of 60° and 120°, respectively. At these azimuthal 
positions, the chordwise flow components are equal, but the radial flow components 
point outward and inward, respectively. The advance ratio was 0.55 (U^ =110 m/sec 
and ftR = 200 m/sec) for both calculations. Figure 2 compares the calculated pres- 
sure distributions at span stations r/R = 0.85, 0.9, and 0.95. The radial flow 
effect is more prominent near the blade tip, which further confirms that three- 
dimensional modeling must be used to calculate rotor-flow phenomena. It is also worth 
noting that for a blade section at the outboard portion of the blade, the negative 
sweep (iJj = 120°) tends to attenuate the expansion at the leading edge and, hence, to 
reduce the supersonic zone and the intensity of the recompression shock. 


Unsteady Effect 

For the same straight-tip rotor blade, figure 3 compares unsteady, quasi-steady, 
and ONERA-measured surface-pressure distributions at the three different span stations 
(r/R = 0.85, 0.9, and 0.95) for azimuthal angles from 30° to 180° at 30° increments. 
The advance ratio was 0.55 (U^ =110 m/sec and QR = 200 m/sec). The agreement 
between unsteady calculations and test data is generally very good. A time retarda- 
tion of the shocks owing to unsteadiness is indicated in the following comparison: 

(1) the shock at 120° azimuth is stronger than that at 90° azimuth in both the 
unsteady calculation and the ONERA measurement; (2) the quasi-steady calculation pre- 
dicts no shock at 150° azimuth, whereas the unsteady calculation catches the shock 
measured in the experiment; and (3) the quasi-steady calculations predict stronger 
shock waves in the first quadrant and weaker shockwaves in the second quadrant. This 
comparison shows that the quasi-steady theory does not adequately predict the strong 
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shock waves that persist over a range of azimuthal angles at extremely high advance 
ratios. Further comparison was made at a lower advance ratio (0.4) for the same 
rotor blade. Figure 4 shows the comparison among unsteady, quasi-steady, and test 
surface pressure distributions. Agreement is quite good for this case in which the 
flows are either subsonic or subcritical with very weak shock waves. It is also 
noted that the quasi-steady theory gives very good pressure predictions for a rotor 
blade near 90° azimuth for flow with moderate shocks and is useful for design work. 


Planform Geometry 

To illustrate the capability of predicting the flow on realistic rotor blades, 
an unsteady calculation was performed for a nonlifting ONERA swept-tip blade 
(ref. 12) at an advance ratio of 0.5 (U^ = 105 m/sec, QR = 210 m/sec). This swept- 
tip blade geometry has a kink at r/R = 0.85, with a 30° sweep angle. Figure 5 shows 
the comparison between unsteady results and ONERA test data. Again, the correlation 
is quite good. The code accounts for the non-straight leading edge and its effect on 
the pressure distributions. Figure 6 shows the evolution with azimuth angle of the 
pressure coefficients at some chordwise points in the same span section. The com- 
parison of full-potential computed results (at x/c =0.3 and r/R = 0.915) and 
ONERA test data (x/c = 0.293 and r/R = 0.911) reveals that full-potential theory, 
even with equations in nonconservative form, is able to predict the transient shock 
position and hence the shock moving speed. However, the shock strength is slightly 
overpredicted, a common result for an inviscid, full-potential calculation. 


CONCLUSIONS 


A finite-difference method for predicting the three-dimensional, unsteady tran- 
sonic flow over a nonlifting helicopter rotor blade in forward flight was presented. 
The method solves the unsteady, transonic, full-potential equations in nonconserva- 
tive form in a blade-fixed mesh system in which the blade surface is conformal with 
a coordinate surface; the system is particularly suitable for modeling the geometry 
effect of the blade. The code based on this method is shown to be a useful tool for 
analyzing unsteady aerodynamics of helicopter rotors. Like its quasi-steady com- 
panions (ref. 11), the unsteady code has several desirable features; for example, 

1. The capability of treating blades of nearly arbitrary geometry 

2. The capability of predicting flow over a blade at any azimuthal angle 

3. The option to restrict calculations to the flow over the outboard portion of 
the blade for computational efficiency 

Computed results obtained from this new code were compared with ONERA data for 
both straight- and swept-tip blades. Good comparison between experimental and 
theoretical pressure distributions of blades could be obtained in transonic flows. 

It is concluded that 

1. The transonic phenomena taking place on the tip of a rotor blade in nonlift- 
ing regime are basically three-dimensional and unsteady 
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2. The quasi-steady theory predicts good unsteady pressure distributions for a 
rotor blade in flow which is either entirely subsonic or subcritical with very weak 
shocks 


3. The quasi-steady theory gives good pressure distributions for a rotor blade 
near 90° azimuth for flow with moderate shock waves; thus, it is useful for design 
work 


4. The unsteady theory is necessary to predict unsteady transonic load for a 
rotor blade in flow that is either supercritical or transonic with moderate shock 
waves 


5- The full-potential equations do not necessarily have to be cast in conserva- 
tion form to catch the transient shock wave movements 

The extension of the method to the lifting rotor blade depends on a good wake 
modeling. For a prescribed wake, the simple model used in reference 13 should be a 
good candidate. For a free wake, a suitable model for transonic potential flow 
equations appears to be a very difficult problem. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, June 13, 1984 
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Figure 1 


{j = 0.335, U^. = 73.7 m/sec, 12 R = 220 m/sec, FLIGHT ALTITUDE = 2000 m 


CONSTANT-MACH LINES 


CONSTANT-SWEEP ANGLE (<;>) LINES 



r- 

CONSTANT-INCIDENCE LINES 



- Constant-Mach number, constant-sweep angle, and constant incidence contours 

on a rotor disk. 
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/j = 0.55, =110 m/sec, flH = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE \ji = 30° 



(a) Azimuth angle 30°. 

fj = 0.55, U^ = 110 m/sec, SIR = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE = 60° 



Figure 3.- Unsteady effect. Comparison among unsteady and quasi-steady calculated, 
and measured, surface-pressure distributions at three span stations near the tip 
of a model rotor blade. 
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fj = 0.55, =110 m/sec, ftR = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE tfr - 120° 


12 T STATION = 1 



STATION - 2 




CENTER OF 
REVOLUTION 


(d) Azimuth angle 120°. 


Figure 3.- Continued. 



H = 0.55, — 110 m/sec, S2R = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE 4/ • 150° 





1 2 3 



(e) Azimuth angle 150°. 

y = 0.55, = 1 10 m/sec, HR = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE * = 180° 






18 


(f) Azimuth angle 180°. 
Figure 3.- Concluded, 



V = 0.4, = 80 m/sec, HH - 200 m/sec 

UNSTEADY 

QUASI STEADY 

• • • ONERA 
AZIMUTH ANGLE v = 60° 



CENTER OF 
REVOLUTION 


(a) Azimuth angle 60°. 

V = 0.4, U^ = 80 m/sec, aR = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE $ = 90° 



(b) Azimuth angle 90°. 

Figure 4.- Assessment of the quasi-steady theory. Comparison among quasi-steady and 
unsteady calculated, and measured, surface-pressure distributions at three span 
stations near the tip of a model rotor blade. 
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fj = 0.4, = 80 m/sec, HR = 200 m/sec 

UNSTEADY 

QUASI-STEADY 

• • • ONERA 
AZIMUTH ANGLE 0 = 120° 



(c) Azimuth angle 120°. 


p = 0.4, U^ = 80 m/sec, S2R « 200 m/sec 

UNSTEADY 

QUASI STEADY 

• • • ONERA 
AZIMUTH ANGLE 0 * 150° 



(d) Azimuth angle 150°. 


Figure 4.- Concluded 



n - 0.5, U oo =105 m/sec, £2R = 210 m/sec 

UNSTEADY 

• • • ONERA 
AZIMUTH ANGLE 0-60° 



(a) Azimuth angle 60°. 

^ = 0.5, U oc = 105 m/sec, HR = 210 m/sec 

UNSTEADY 

• • • ONERA 
AZIMUTH ANGLE ^ = 90° 



(b) Azimuth angle 90 °. 

Figure 5.- Planform effect. Comparison between unsteady calculated and measured, 
surface-pressure distributions at three span stations near the tip of a swept-tip 
(back 30°) model rotor blade. 



fj = 0.5, = 105 m/sec, TiR = 210 m/sec 

UNSTEADY 

• • • ONERA 
AZIMUTH ANGLE $ = 120° 





CENTER OF 
REVOLUTION 



(c) Azimuth angle 120°. 


jj = 0.5, U^ = 105 m/sec, HR = 210 m/sec 

UNSTEADY 

• • • ONERA 
AZIMUTH ANGLE ^ = 150° 





CENTER OF 
REVOLUTION 



(d) Azimuth angle 150°. 


Figure 5.- Continued 


*i = 0.5, = 105 m/sec, 12 R = 210 m/sec 

UNSTEADY 

• • • ONERA 
AZIMUTH ANGLE v = 180° 





CENTER OF 
REVOLUTION 



(e) Azimuth angle 180°. 
Figure 5.- Concluded. 


jj = 0.5, U^ = 105 m/sec, OR = 210 m/sec 

PRESENT METHOD 

ONERA TEST DATA 







PRESSURE SENSORS ON THE UPPER SURFACE 
OF THE ROTOR BLADE SECTION AT r/R = 0.9147 


Figure 6.- Azimuthal pressure evolution on a swept-tip rotor blade 
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